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Abstract 
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Because glueballs are SU(3)pi avor singlets, they are expected to couple equally to 

u, d, and s quarks, so that equal coupling strengths to ir + -K~ and K + K~ are predicted. 

However, we show that chiral symmetry implies the scalar glueball amplitude for Go — ► 

, qq is proportional to the quark mass, so that mixing with ~ss mesons is enhanced 

Oh, and decays to K + K~ are favored over ir + ir~. Together with evidence from lattice 

Qh. calculations and from experiment, this supports the hypothesis that /o(1710) is the 

^ . ground state scalar glueball. 
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1. Introduction. The existence of gluonic states is a quintessential prediction of 
Quantum Chromodynamics (QCD). The key difference between Quantum Electrodynamics 
(QED) and QCD is that gluons carry color charge while photons are electrically neutral. 
Gluon pairs can then form color singlet hadronic bound states, "glueballs," like mesons and 
baryons, which are color singlet bound states of valence quarks. jT Because of formidable 
experimental and theoretical difficulties, it is frustrating, though not surprising, that this 
simple, dramatic prediction has resisted experimental verification for more than two decades. 
Quenched lattice simulations predict that the mass of the lightest glueball, G , a scalar, is 
near ~ 1.65 GeV,[2] but the prediction is complicated by mixing with qq mesons that require 
more powerful computations. Experimentally the outstanding difficulty is that glueballs are 
not easily distinguished from ordinary qq mesons, themselves imperfectly understood. This 
difficulty is also exacerbated if mixing is appreciable. 

The most robust identification criterion, necessary but not sufficient, is that glueballs 
are extra states, beyond those of the qq meson spectrum. This is difficult to apply in practice, 
though ultimately essential. In addition, glueballs are expected to be copiously produced in 
gluon rich channels such as radiative J / ip decay, and to have small two photon decay widths. 
These two expectations are encapsulated in the quantitative measure "stickiness," [3| which 
characterizes the relative strength of gluonic versus photonic couplings. 

Another popular criterion is based on the fact that glueballs are S77(3)Fi av or singlets 
which should then couple equally to different flavors of quarks. However we show here 
that the amplitude for the decay of the ground state scalar glueball to quark-antiquark is 
proportional to the quark mass, /A (Go — > qq) oc m q , so that decays to ss pairs are greatly 
enhanced over uu + dd, and mixing with ss mesons is enhanced relative to uu + dd. We 
exhibit the result at leading order and show that it holds to all orders in standard QCD 
perturbation theory. 

The result has a simple nonperturbative physical explanation, similar, though different 
in detail, to the well known enhancement of tt — > \w relative to ir — > ev. For m q = chiral 
symmetry requires the final q and q to have equal chirality, hence unequal helicity, so that 
in the Go rest frame with z axis in the quark direction of motion, the total z component of 
spin is nonvanishing, \Sz\ = 1. Because the ground state Go gg wave function is isotropic 
(L = S = 0), the qq final state is pure s- wave, 2 L = 0. The total angular momentum is 
zero, and since there is no way to cancel the nonvanishing spin contribution, the amplitude 
vanishes. With one power of m q ^ 0, the q and q have unequal chirality and the amplitude 
is allowed. 

The enhancement is substantial, since m s is an order of magnitude larger than m u and 

integrating over the gluon direction to project out the s-wave gg wave function is equivalent to integrating 
over the final quark direction with the initial gluon direction fixed. 
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But for scalar glueballs of mass ~ 1.5 — 2 GeV, r(G — > ss), is suppressed of order 
(m s /mG ) 2 , so that it may be smaller than the nominally higher order Go — > qqg process, 
which is SU(3)fiscvot symmetric. We find that the soft and collinear quark-gluon singularities 
of Go — > qqg vanish for m q = 0, as they must if Go — > qq is to vanish at one loop order for 
m q = 0. Unsuppressed, flavor-symmetric Go — > qqg decays are dominated by configurations 
in which the gluon is well separated from the quarks, which hadronize predominantly to 
multi-body final states. Go may also decay flavor symmetrically to multigluon final states 
(n > 3), via the three and four gluon components of G 2 in equation (1), by higher dimension 
operators that arise nonperturbatively, or by higher orders in perturbation theory. The IR 
singularities of Go — > ggg are cancelled by virtual corrections to the Go wave function, while 
configurations with three (or more) well separated gluons hadronize to multihadron final 
states. The enhancement of ss relative to uu + dd is then most strongly reflected in two 
body decays: we expect K + K~ to be enhanced relative to 7r + 7r~, while multibody decays 
are more nearly flavor symmetric. 

Glueball decay to light quarks and/or gluons cannot be computed reliably in any fixed 
order of perturbation theory. However, the predicted ratio, T(Go — > ss)/T(Go — > uu + dd) 
1, is credible since it follows from an analysis to all orders in perturbation theory and, in 
addition, from a physical argument that does not depend on perturbation theory. The 
implication that T(Go — > K + K~) T(Go — > 7r + 7r~) is less secure and is best studied on the 
lattice. Remarkably, it is supported by an early quenched study of G decay to pseudoscalar 
meson pairs for two "relatively heavy" S{/(3)Fiavor symmetric values of m q , corresponding 
to mps — 400 and ~ 630 MeV.jHj Linear dependence on m q implies quadratic dependence 
on mps,jE] which is consistent at la with the lattice computations. [5] Chiral suppression 
could then be the physical basis for the unexpected and unexplained lattice result. With 
subsequent computational and theoretical advances in lattice QCD, it should be possible 
today to verify the earlier study and to extend it to smaller values of mps, nearer to the 
chiral limit and to the physical pion mass. If the explanation is indeed chiral suppression, 
then the couplings of higher spin glueballs should be approximately flavor symmetric and 
independent of mps, a prediction which can also be tested on the lattice. 

Enhanced strange quark decay changes the expected experimental signature and sup- 
ports the hypothesis that / (1710) is predominantly the ground state scalar glueball. This 
identification was advocated by Sexton, Vaccarino, and Weingarten, [5] and is even more 
compelling today in view of recent results from J /if) decay obtained by BES[7J|H] - see P 
for an overview of the experimental situation. 

In section 2 we compute *M(Go — > qq) at leading order for massive quarks, with the 
expected linear dependence on m q . In section 3 we show that A^(Go — > qq) oc m q to any 
order in as- In section 4 we describe the infrared singularities of M.{Gq — > qqg). We conclude 
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with a brief discussion, including experimental implications. 

2. G — > qq at leading order. Consider the decay of a scalar glueball G with 
mass Mg to a qq pair 3 with quark mass m q . The effective glueball-gluon-gluon coupling is 
parameterized by 

£eff = foG G aflu Ga U (1) 

where Go is an interpolating field for the glueball, G aiiu is the gluon field strength tensor with 
color index a, and f is an effective coupling constant with dimension 1/M that depends on 
the Gq wave function. The gg — > gg scattering amplitude can be written as 

M(gg -> gg) = e^M^igg -> gg) (2) 

where ej M = e^(pi, Xi) with i — 1,2 are the polarization vectors for massless constituent 
gluons with four momentum and polarization Aj. Using equations (1) and (2), summing 
over the polarizations Aj, and averaging over the gluon direction in the G rest frame to 
project out the s-wave, we obtain 

M(G -gg) = ^/ cMX^M^gg - gg), (3) 

where = 2p^ - M G g» v projects out the |(++) + ( ) > helicity state that couples 

to Ga^G^ in equation (1). 

From the lowest order Feynman diagrams we obtain 

X^M, V = J^J^L _ _ s(ps> ^ Mp4 , ^ (4) 

where W3, t>4 are the q, q spinors for quark and antiquark with center of mass momenta P3,Pa, 
helicities h^, h±, color indices i,j and center of mass velocity j3. Equation (4) includes a color 
factor from the color singlet gg wave function, 

" ^TT " ~Y ij [b) 

Performing the angular integration, the decay amplitude is 

KAin — \ f 1QkV2 m q l + /3_ 

M{G -> gg) = -foa s — — lo S]— g u 3 v 4 d i:j . (6) 



3 Elastic scattering, gg — > gg, contributes to the glueball wave function, not to the decay amplitude. 
The dissociation process, gg — > qq + qq, in which each gluon makes a transistion to a color-octet qq pair, 
is kinematically forbidden for m q ^ 0; an additional gluon exchange is required to allow it to proceed 
on-mass-shell, which is therefore of order g 4 in the amplitude. 
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Squaring equation (6), summing over quark helicities and color indices, and performing the 
phase space integration, the decay width is 

T(Go - qq) = ijp o? s flm\M G (5 log 2 [±|. (7) 

Notice that an explicit factor m q appears in the gg —>■ qq amplitude, equation (4), 
which is not averaged over the initial gluon direction and which clearly has contributions 
from higher partial waves, J > 0. It may then appear that chiral suppression applies not 
just to spin glueballs but also to glueballs of higher spin. However when equation (4) is 
squared and the phase space integration is performed, a factor 1/m 2 , results from the t and u 
channel poles, which cancels the explicit factor m q in the numerator, so that the total cross 
section a(gg — > qq) does not vanish in the chiral limit, because of the J > partial waves. 

3. Go — > qq to all orders. - We now show that M.{G$ — > qq) vanishes to all orders in 
perturbation theory for m q = 0. Consider the Lorentz invariant amplitude 

M x (p 1 ,p 2 ,p 3lP4 )=X^M llu (8) 

where M.^ is defined in equation (2) and X^ LU below (3). The perturbative expansion for 
M.x is a sum of terms arising from Feynman diagrams with arbitrary numbers of loops. 
After evaluation of the loop integrals, regularized as necesary, Adx is a sum of terms, 

M x = J2 1I (p3,X3)Tiu(p 4 ,X'i), (9) 

i 

where Us, are respectively massless fermion and antifermion spinors[TU] of chirality X3 ; X4- 
The Ti are 4x4 matrices, each a product of momentum-contracted Dirac matrices, 

Ti = fcil fci2 ■ ■ ■ fci ni (10) 

where each ki a is one of the external four-momenta, Pi,P2,P3,Pa- 

Chiral invariance for m q = implies that the number of factors, rii, in equation (10) is 
always odd. Since all external momenta vanish and the spinors obey /P3U3 = /p^Ui = 0, by 
suitably anticommuting the faa, each term in equation (9) can be reduced to a sum of terms 
linear in j6i and j^ 2? which we choose to be symmetric and antisymmetric, 

U(P3, X3)IYu(p 4 , Xi) =u(P3,X3)[Si(s,t,u)(ft 1 + fa) + A^S, t, U) ( fry- ^ 2 )]«(P4, Xi) • ( U ) 

The coefficients A iy Si are Lorentz invariant functions of the Mandelstam variables s,t,u. 
Since pi + P2 = P3 + Pi, the symmetric term vanishes and equation (9) reduces to 

M x = A(s,t,u) u(p3,Xs)(fti- h)u(P^Xi)- (12) 
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where A(s, t, u) = J2i Ai(s, t, u). 

Next consider the integration over the gluon direction, equation (3). In the G rest 
frame with the z-axis chosen along the quark direction of motion, z = j3 3 , we integrate 
over dfl = d 2 p\ = d(j)idcos6\. The Mandelstam variables are then s = Mq and u,t = 
— ^Mq(1 ± cos^i). Since the color and helicity components of the Go wave function are 
symmetric under interchange of the two gluons, Bose symmetry requires A(s, t, u) to be odd 
under pi <-> p 2 . In our chosen coordinate system A is a function only of cos#i, and must 
therefore be odd, A(— cos^) = — ^(cos^!). But evaluating ^3(^1— /p<i)u± explicitlv[TU] we 
find 

h>A = M 2 G e-** 1 sin^ (13) 

which is even in cos#i, while the azimuthal factor, e _l<?il , provides the required oddness under 
Pi *-* P2- e - *^ 1 — * e~ l ^ 1+7T ^ = — e - *^ 1 . Consequently the integral / dcosOiA vanishes, and 
Ai(Go — > qq) = to all orders in the chiral limit. In fact, because of our choice of axis, 
z = P3, the integral over <pi also vanishes. For other choices of z the azimuthal and polar 
integrals do not vanish separately, but the full angular integral, / d 2 p\, vanishes in any case. 

For nonvanishing quark mass, m q 7^ 0, chirality-flip amplitudes contribute. With one 
factor of m q from the fermion line connecting the external quark and antiquark, the Tj 
matrices in equation (9) include products of even numbers of Dirac matrices, i.e., in 
equation (10) may be even. Beginning in order m q there are then nonvanishing contributions 
to Ai(Go — > qq), like the leading order term shown explicitly in equation (6). 

The vanishing azimuthal integration for z = p3 reflects the physical argument given in 
the introduction. The factor e _l<?il corresponds to Sz = 1 from the aligned spins of the q and 
q, while the absence of a compensating factor in A is due to the projection of the orbital 
s-wave by the / d 2 p\ integration and the absence of spin-polarization in the initial state. 

4. Infrared singularities of Go — > qqg. Although it is of order a|, r(G — ► qqg) 
is not chirally suppressed and may therefore be larger than r(G — > ss), which is of order 
a 2 s x m 2 s /m 2 G . Setting m q = we evaluated the 13 Feynman diagrams using the helicity 
spinor method [.G]j with numerical evaluation of the 9 dimensional integral: 

T(G ^qqg) = Eh 3 ,fc4,A 5 /ps \M(G -^<? 3 ?4#5) I 2 

= ifs T. hiMM Jps / dSl x f dn[ e^X^M^igm ^ q 4 q^) 

x e h pX' ap M^{g' ig ' 2 ^ ? A q 3 g 5 )*. (14) 

Details will be presented elsewhere. jTTj We focus here on the infrared singularities, which 
provide a consistency check at one loop order that Ai(Go — > qq) vanishes in the chiral limit. 

In general there could be soft IR divergences for E q , E q , E g — > and collinear divergences 
for 9 qg , 9 qg , 9 qq — > 0. In fact, only the qq collinear divergence occurs, as can be seen from the 
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Figure 1: Distributions for Go — > in arbitrary units. In figure (la) the dot-dashed line 
is dN/dEg and the dashed line is dN/dE q . In figure (lb) the dot-dashed line is dN/cos9 qg 
and the dashed line is dN/dcos9 qq . 

distributions in figure 1, obtained by imposing only the cut 9 qq > 0.1 in the Go rest frame: 
neither dN/dE g nor dN/dE q diverge at low energy, and only dN/9 qq diverges at 9 qq — > 0. 
Instead dN/dE g diverges at the maximum energy, E g = mc /2, and dN/d9 qg diverges for 
9 qg — > it. Both of these divergences are kinematical reflections of the collinear singularity at 
9q q — > 0, for which the qq pair with m qq = recoils with half of the available energy against 
the gluon in the opposite hemisphere. 

This is precisely the pattern of divergences required if Go — > qq is chirally suppressed to 
all orders and, in particular, at one loop. For if there were soft divergences in any of E q , Eg, E g 
or collinear divergences in 9 qg and 9 qg , then the resulting singularities at m qg , m qg — > would 
have to be cancelled by virtual corrections to Go — > qq, such as gluon self energy contributions 
to the quark propagator. The absence of these singularities is a consistency check (i.e., a 
necessary condition) that Go — > qq is chirally suppressed at one loop order. The collinear 
divergence for 9 qq — > is cancelled by quark loop contributions to the gg — > gg amplitude, 
which in the present context are one loop corrections to the Go wave function. 

5. Discussion. We have shown to all orders in perturbation theory and with a 
simple, nonperturbative physical argument that the ground state J = glueball has a 
chirally suppressed coupling to light quarks, -M(G — > qq) oc m q , with corrections of higher 
order in m q /rriG- From equation (7) with m u ,md, m s varied within la limits, jlj T(Go — > ss) 
dominates T(Go — > uu + dd) by a factor between 20 and 100. Flavor symmetry is reinstated 
for Go — > qqg when the gluon is well separated from the q and q. For sufficiently heavy nic 
one can test this picture by measuring strangeness yield as a function of thrust or sphericity, 
with enhanced strangeness in high thrust or low sphericity events, but it is unclear if this 
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is feasible for mg ~ 1.7 GeV. It is more feasible for the ground state pseudoscalar glueball, 
which is expected to be heavier than the scalar and which we also expect to be subject to 
chiral suppression. 

For light scalar glueballs, the best hope to see strangeness enhancement is the two body 
decays, G — > K + K~ . Since G — > qqg and Go —* 999 are not chirally suppressed, 
naive power counting suggests r(G — > qqg + ggg) > r(G — > ss), so that n > 3 parton decay 
amplitudes are probably the dominant mechanism for multihadron production. Then KK 
will dominate two body decays while multiparticle final states are approximately SU( 3) Flavor 
symmetric, up to phase space corrections favoring nonstrange final states. 

Chiral suppression has a major impact on the experimental search for the ground state 
scalar glueball. Candidates cannot be ruled out because they decay preferentially to strange 
final states, especially KK, and mixing with ss mesons may be enhanced. This picture of 
a chirally suppressed Go fits nicely with the known properties of the /o(1710) meson. It is 
copiously produced in radiative ip decay in the ip —>■ ^KK channel[7j and in the gluon-rich 
central rapidity region in pp scattering. |12j has a small 77 coupling,^] has a mass consistent 
with the prediction of quenched lattice QCD, 2J and has a strong preference to decay to 
KK, with B(nn)/B(KK) < 0.11 at 95% CL. 8] As a rough estimate of the stickiness, [3] we 
combine the 77 95% CL upper limit with central values for ip radiative decay [Zj, with the 
result S(/ (1710)) : 5(^(1525)) : 5(/ 2 (1270)) ~ (> 36) : 12 : 1. A more complete discussion 
of the experimental situation will be given elsewhere |llj — see also jHj- 

The interpretation of / (1710) as the chirally suppressed scalar glueball can be tested 
both theoretically and experimentally. Lattice QCD can test the prediction that Go — > KK is 
enhanced for the ground state J = glueball but not for J > 0. With an order of magnitude 
more J/ip decays than BES II, experiments at BES III and CESR-C will extend partial 
wave analysis to rarer two body decays and to multiparticle decays. They could confirm 
B(fo — > KK)/B(f — > 7T7r) ^> 1 and, if, as is likely, the rate for multiparticle decays is big, 
the lower bound on B{ip — > 7/o(1710) will increase beyond its already appreciable value from 
KK alone. A large inclusive rate B(ip — > 7/0), a large ratio B(f — ► KK)/B(f — > nn), 
and approximately flavor symmetric couplings to multiparticle final states would support 
the identification of /o(1710) as the chirally suppressed, ground state scalar glueball. 

Note added: Cornwall and Soni observed that the trace anomaly [15^ implies chiral sup- 
pression of the Go coupling to tttt and KK at zero four-momentum; however, the extrapo- 
lation to the G mass-shell, e.g., at ~ 1700 MeV, is not under control and could be large. 
Similarly, in the conjectured AdS/CFT approach to QCD the scalar glueball is a dilaton. [To] 
which therefore has chirally suppressed two body couplings at zero momentum. I thank A. 
Soni, H. Murayama, and M. Schwartz for bringing this work to my attention. 
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